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Abstract 



■ We consider cosmological solutions and their stability with respect to homoge- 

neous and isotropic perturbations in the braneworld model with the scalar-curvature 
term in the action for the brane. Part of the results are similar to those obtained 
by Campos and Sopuerta for the Randall-Sundrum braneworld model. Specifically, 
| the expanding de Sitter solution is an attractor, while the expanding Friedmann 

solution is a repeller. In the braneworld theory with the scalar-curvature term in 
the action for the brane, static solutions with matter satisfying the strong energy 
condition exist not only with closed spatial geometry but also with open and flat 
ones even in the case where the dark-radiation contribution is absent. In a certain 
^ ' range of parameters, static solutions are stable with respect to homogeneous and 

bJQ' isotropic perturbations. 



PACS number(s): 04.50. +h 

1 Introduction 



Braneworld models have become quite popular in the high-energy and gravitational phy- 
sics during last several years. There are at least three reasons for this. Firstly, most of 
the brane models, in particular, the Randall-Sundrum model |T], are inspired directly by 
the popular string theory (specifically, by the Hofava-Witten model j2]) and represent an 
interesting alternative for compactification of extra spatial dimensions. Secondly, there 
remains a possibility of solving the mass-hierarchy problem in the context of large extra 
dimensions (see |H1 E]). And thirdly, it has been shown that models with non-compact 
extra dimensions and branes can be consistent with the current gravity experiments. 

Braneworld models also turned out to be consistent with modern cosmology, at the 
same time exhibiting new specific features. First of all, this concerns the braneworld 
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models of "dark energy," or of the currently observed acceleration of the universe (see 
jl] in this respect), but this is also true with respect to the early stages of cosmological 
evolution such as inflation (see, e.g., jH]). 

A broad class of cosmological solutions in braneworld theory was systematically an- 
alyzed by Campos and Sopuerta [B] using convenient phase space variables similar to 
those introduced in These authors gave a complete description of stationary points 
in an appropriately chosen phase space of the cosmological setup and investigated their 
stability with respect to homogeneous and isotropic perturbations. The authors worked 
in the frames of the Randall-Sundrum braneworld theory without the scalar-curvature 
term in the action for the brane. In this paper, we are going to extend this analysis by 
investigating the case where this term is present in the action. 

The scalar-curvature term for the brane was introduced in JSJ (see also ||) as a method 
of making gravity on the brane effectively four- dimensional even in the flat infinite bulk 
space. The corresponding cosmological models were initiated in 0EH]. The necessity of 
this term on the brane arises when one considers the complete effective action of the theory 
(see, e.g., [II]). From this viewpoint, the term with brane tension represents the zero-order 
term in the expansion of the total action for the brane in powers of curvature of the induced 
metric, while the next term in the action for the brane is exactly the scalar-curvature 
Hilbert-Einstein term. From another (complementary) viewpoint jH], this contribution is 
induced as a quantum correction to the effective action for the brane gravity after one 
takes into account the quantum character of the matter confined to the brane. 

In this paper, we will show that the presence of the scalar-curvature term in the action 
for the brane enlarges the class of cosmological solutions, in particular, it allows for the 
existence of static cosmological solutions with ordinary matter content even in the case 
of spatially flat or open cosmology, a specific feature of the braneworld cosmology that 
was also noted in jHJ H2 H] • In the Randall-Sundrum model, flat or open static universe 
with ordinary matter requires the presence of the negative dark-radiation term, while, in 
the model with scalar-curvature term in the action for the brane to be discussed in this 
paper, the presence of dark radiation may be not necessary. 

The paper is organized as follows: In Sec. we introduce the basic equations de- 
scribing the theory with one brane in the five-dimensional bulk space. Then, in Sec. 
we introduce a method for the investigation of static solutions. Using this method, we 
describe both the general-relativity case and the case of braneworld theory that contains 
the scalar-curvature term in the action for the brane. In the partial case of spatially flat 
universe, we determine the range of parameters where static solutions exist and where 
they are stable with respect to homogeneous and isotropic perturbations. In Sec. El we 
develop an approach which is used to describe non-static solutions, and apply it to the 
braneworld theory under investigation. In Sec. we return to the issue of static solutions 
with arbitrary spatial curvature using conveniently chosen variables. Finally, in Sec. HI we 
present our conclusions. 
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2 Basic equations 



In this paper, we consider the case where the braneworld is the time-like boundary of a 
five-dimensional purely gravitational Lorentzian space (bulk), which is equivalent to the 
case of a brane embedded in the bulk with Z2 symmetry of reflection with respect to the 
brane. The theory is described by the action j^j: 



S 



(ft-2A b )-2 / K 

bulk J brane 



(m 2 R - 2a) + 



brane 



L(h. 



ab ■■ 



(2.1) 



brane 



Here, 1Z is the scalar curvature of the metric g ab in the five-dimensional bulk, and R is 
the scalar curvature of the induced metric h ab = g ab — n a n b on the brane, where n a is the 
vector field of the inner unit normal to the brane, and the notation and conventions of 
[T3~] are used. The quantity K = K ab h ab is the trace of the symmetric tensor of extrinsic 
curvature K ab = h^V c n b on the brane. The symbol L(h ab ,(p) denotes the Lagrangian 
density of the four-dimensional matter fields whose dynamics is restricted to the brane 
so that they interact only with the induced metric h ab . All integrations over the bulk and 
brane are taken with the natural volume elements ^—gd 5 x and \J — hd 4 x, respectively, 
where g and h are the determinants of the matrices of components of the corresponding 
metrics in a coordinate basis. The symbols M and m denote, respectively, the five- and 
four-dimensional Planck masses, Ab is the bulk cosmo logical constant, and a is the brane 
tension. 

The term containing the scalar curvature of the induced metric on the brane with the 
coupling m? in action ()2.1|) was originally absent from braneworld models. However, very 
soon it became clear that this term is qualitatively essential for describing the braneworld 
dynamics and, moreover, it is inevitably generated as a quantum correction to the matter 
action in ()2.1|) — in the spirit of an idea that goes back to Sakharov [T^j (see also [T5j). 
Note that the effective action for the brane typically involves an infinite number of terms 
of higher order in curvature (this was pointed out in [Hj for the case of braneworld theory; 
a similar situation in the context of the AdS/CFT correspondence is described in In 
this paper, we retain only the terms linear in curvature. The effects of the curvature term 
on the brane in linear approximation were studied in d] , where it was shown that it 
leads to four- dimensional law of gravity on sufficiently small scales. The presence of this 
term also leads to qualitatively new features in cosmology, as demonstrated, in particular, 
in E] • For some recent reviews of the results connected with the induced-gravity term 
in the brane action, one may look into [17 . 

Variation of action (j2.1j) gives us the equations describing the dynamics in the bulk 



Gab + Ab^ab = 



(2.2) 



and on the brane 



m 2 G ab + ah ab = r ab + M 3 (K ab - Kh a 



(2.3) 
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The second equation generalizes the Israel junction conditions and the corresponding 
equation of the Randall-Sundrum model to the presence of the scalar-curvature term on 
the brane (m^O). 

One can consider the Gauss identity 

Rab/ = hJh h a h c k h d -R hk i + K ac K b d - K bc K a d , (2.4) 

and, following the original procedure of contract it once on the brane using equations 
(|2.2|) and (|2.3jl . In this way, one obtains the effective equation on the brane that generalizes 
the result of to the presence of the brane curvature term: 

G ab + A cS h ab = ^r ab + 1 + 2 ^ 2/3M6 (M*Q ab - W ab ) , (2.5) 

where 

Ars . A b o 2 2 o 3M 6 

Acff = 1 + 2amV3M« ' RS = T UP ' Mp = m + ~2a~ ' (2 ' 6) 

Qafe = 7^EE ab — E ac E c b + - ^E cd E cd — t^E 2 ^ h ab (2.7) 

is the quadratic expression with respect to the 'bare' Einstein equation E ab = m 2 G ab — r ab 
on the brane, E = h ab E ab , and W ab = h c a h e b W c d e fn d n^ is a projection of the bulk Weyl 
tensor W abc d to the brane. The limit of m — > leads us to the original result of [T§] . 

Contracting Eq. (|2.5j) once again, one obtains the following scalar equation on the 
brane: 

M\R - 2A b ) + (m 2 G ab + ah ab - r ab ) (m 2 G ab + ah ab - r ab ) - \{m 2 R - 4a + rf = , (2.8) 

where r = r ab h ab . It is argued in that, if one is interested only in the evolution on 
the brane and if no additional boundary or regularity conditions are specified in the bulk 
(as often is the case), then one is left only with this scalar equation. Having obtained 
a solution of (|2.8|) . one can integrate the gravitational equations in the bulk using, for 
example, Gaussian normal coordinates [THj I2()j. 

The brane filled by a homogeneous and isotropic ideal fluid is described by the Fried- 
mann-Robertson-Walker (FRW) metric 



ds 2 = -dt 2 + a 2 (t) 



dr2 +r 2 dn 2 



(2.9) 



1 — nr 2 

where k = —1, 0, or +1. The stress-energy tensor of the ideal fluid takes the form 

Tab = (P + P)UaUb + Phab , (2.10) 

where p is the energy density, and p is the pressure. In this paper, we assume a simple 
equation of state p = wp, where w is constant. 
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Equation ()2.8|) leads to the following cosmological equation in the braneworld theory: 



9M e 



° +H 2 + k _K 



+ 



p(l + 3w) -2cr + 6m 2 - 

a 



p + a -3m z [H z + — 



(2.11) 

where if = a/a is the Hubble parameter. After multiplying this equation by a 3 d and using 
the conservation of energy 

p + 3H(p + p) = 0, (2.12) 



one can integrate Eq. (|2.11J) and obtain the first-order differential equation 2 

2 



I i 7-7-2 K p + O" 

m 4 ( if 2 + — 

a z 



3m 2 



K A>- 



C 



a z a 4 



(2.13) 



where C is the integration constant corresponding to the black-hole mass of the Schwarz- 
schild-(anti)-de Sitter solution in the bulk and associated with what is called "dark radi- 
ation." 

Equation (J2.13|) can be solved with respect to the Hubble parameter: 



H 2 + 



K 



p + a + 2 



3m 2 



i 2 



i± 4/1 + 



p 



a 



3m 2 



Ab 
6 



C 



(2.14) 



where i = 2m 2 /M 3 is the crossover length scale characterizing the theory. The £ ±' signs 
give the two different branches of braneworld solutions which, to some extent, are con- 
nected with the two different ways in which the Schwarzschild-(anti)-de Sitter bulk space 
is bounded by the brane: the inner normal to the brane can point either in the direction 
of increasing or decreasing of the bulk radial coordinate. Following [4 , we refer to the 
model with the lower sign (' — ') as to BRANE1, and to that with the upper sign ('+') as 
to BRANE2. 

In the formal limit of m — > 0, the BRANE1 model reduces to the well-known expression 
obtained in the Randall-Sundrum case 



H 2 + 



(P + g) 
9M 6 



Ab C 
6 a 4 



(2.15) 



while the BRANE2 model does not exist in this limit. This can be viewed as a simple 
consequence of the fact that Eq. (|2.13|) . being quadratic with respect of the Hubble pa- 
rameter H in the case m^O, reduces to a linear equation with respect to H in the formal 
limit m —>■ 0. 



3 Static solutions 

We proceed to the investigation of static solutions in the braneworld theory. First, we 
review the situation in the general relativity theory, where the static Einstein's universe 
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is known to be unstable. We apply the same method for the investigation of the stability 
of static solutions in the braneworld theory and demonstrate new possibilities arising in 
this theory. Firstly, unlike general relativity, the braneworld theory admits static solutions 
even in the case of spatially flat or open universe, which, in the context of the Randall- 
Sundrum model, was also noted in in the Randall-Sundrum model, static spatially 
flat or open universes with matter satisfying the strong energy condition are possible only 
with negative value of the dark-radiation constant C, while, in the presence of the scalar- 
curvature term in the action for the brane, a spatially flat or open universe can be static 
even when this constant is zero. Secondly, these solutions can be stable in certain region 
of parameters. 



3.1 The case of general relativity 

In this section, we review the situation in the general relativity theory with the standard 
action 

S G r = -ttt; f( R - ^)V^gd 4 x + [ L matterV /^ 4 x (3.1) 



16ttG . 

and describe our method for investigation of the stability of static solutions on the example 
of this theory. A homogenous and isotropic universe with density p and pressure p = wp 
is described by the Friedmann equations 



k 8irGp A 



a 3 3 

giving the following conditions of static solution: 

k SnGpo A 
Oq 3 3 



a AttGo A . 

' 'l + 3w) + ~, (3.3) 



(3.4) 



^(1 + 3.) -1 = 0. (3.5) 
Introducing small homogeneous deviations from stationarity, 

a = ao(l-e), e = e(t) < 1 , (3.6) 

and linearizing the Friedmann equations ()3.2|) and (j3.3j) with respect to the small quantity 
e(t), we have: 

2fte 8tiG(p — p ) 



3 



(3.7) 



We, W= K{1 \ 3W) . (3.* 



6 



We note that the weak energy condition and casuality requirements restrict the value of 
w to lie in the interval w G [—1, +1] (see, e.g., [21]) and, as a result, n = +1. Taking this 
into account, one can see from (|3.8|) that the static solution in this theory is unstable for 
w G [— |, l] , and stable for w G (—1, — ~). 



3.2 The case of braneworld theory 

One can use the method outlined in the previous section to investigate the braneworld 
theory. In this section, we find (analytically and graphically) the regions of existence of 
static solutions in the parameter space and then further restrict the space of parameters by 
the requirement that the described theory be physically compatible with the observations. 



3.2.1 Static solutions and their perturbation 

In the static case, equations ()2.11|) and (|2.13|) give 

9M 6 f^-^\ + [p{\ + 3w) - 2a] [p + a - 3m 2 ^) = (3.9) 

and 

m U-l^J =M {^'Y-^) ' (3 - 10) 

respectively. One should note that Eq. (j3.10|) cannot be obtained by integrating Eq. (|2.11|) 
in the static case because, before integration, Eq. (j2.11|) is multiplied by a, which, in this 
case, is identically zero. Nevertheless, it can be obtained from the original equations 
on the brane (J2.3)) considering the embedding of the brane in the bulk. In fact, this 
equation determines the value of the black-hole constant C for which a solution of (13. 9|) 
is embeddable in the Schwarzschild-(anti)-de Sitter bulk. 

Looking for static solutions and investigating their stability, we proceed as in the pre- 
vious case of general relativity. First, we solve Eq. (j2.11j) with respect to a /a for a = 0: 

q _ 9M 6 (I - f ) + W + 3w) -2a](p + a- 3m 2 ^) 
a 9M 6 + 6m 2 (p + a - 3m 2 ^) 

Studying homogeneous perturbations of this solution, we set a = as the initial condition. 
If a solution is stable (unstable) with respect to such perturbations, it will also be stable 
(unstable) with respect to more generic perturbations with a ^ 0. Thus, we assume that 
the relation between p and a is given by Eq. (j3.10|) at the initial moment of time. 
We express the energy density p through a using Eq. (|3.1(Jj) . 
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and substitute it to ([3.110 : 

a A\ 



a 9M 6 ^18m 2 M 3 A /4-4^-4 

1 V a 6 a 4 



(3.13) 



where 



A 2 = -3a(l + w) + 3(1 + 3u>)m 2 -^ T 3(1 + 3w)M\ ^ - -± - — . (3.15) 

cr V a o a 4 

After that, we introduce the perturbation of the static solution 

a = ao(l-e), e(t)«l, (3.16) 

where a is a stationary solution of ([3.13[l . and expand Eq. ([3.13(1 around the stationary- 
point: 

e = We, (3.17) 



where 



W = Bl (3.18) 

I _|_ 2m? I _k_ _ Aj, _ _C 



i?i = ^ + (l + 3 W ) ft-^r) T5 2 ±5 3 ±54, 



5, 



a \ a a 

crfl + w) /« 2C 



M 3 /4 -4^- 4 V a o a 



6 ^ 



'0 



m 2 2k /k A b C 

2 w» 2C 



(3.19) 



M 3 a 2 , /k _ Ab _ c ' 

6 4 

The lower and upper signs in these equations refer to the BRANE1 and BRANE2 models, 
respectively. 

It can also be noted \A\ that, unlike the general relativity theory, the braneworld theory 
allows for empty (p = p = 0) static solutions. The radius (scale factor) of such a universe 
is given by the expression jl] 

2 k ( 3 am 2 



The spatially flat case requires Ars = 0. 
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3.2.2 Physical restrictions on the constants of the theory 



Observations indicate that, at the present cosmological epoch, our universe expands with 
acceleration, so that the simplest asymptotic condition for such a universe is that the 
energy density will decrease monotonically to zero. Thus, we can consider the low-energy 
behavior of the theory as p — ► 0: 



K 



— = C X + C 2 p. 
a 1 



Ci > , C 2 > . 



which is valid under the condition 

£ 2 p 
3m 2 



< 1 + 



a 



Ab 

3m 2 6 



(3.21] 



(3.22) 



and where we have assumed that the dark-radiation term (with the constant C in (|2.14jl ) 
has negligible contribution. Assuming regime (|3.21j) is realized today, we can write the 
following minimal physical restrictions on the constants of the theory: 



BRANE1: The inequalities 
Ci 



a 2 



3m 



1 + 



a 





3m 2 6 



> 



and 



Co 



1 



3m 2 



1 + ^ 2 fe-f) 



> 



give 



cr > , a > 
BRANE2: The inequalities 



A b m 2 



ArsS T + 3^ >0 



(3.23) 



(3.24) 



(3.25) 



a 2 



3m 2 



1 + A 1+P 



a 



At 

3m 2 6 



> 



and 



C 2 



3m 2 



imply 



a > 



6m 



1 + ^ 2 fe-t)J 
A b m 2 3m 2 



> 



a > 



e 2 



(3.26) 



(3.27) 



(3.28) 



From these relations, one can see that BRANE2 model with the asymptotic conditions 
(|3.21j) . in principle, can have negative brane tension, while the tension of the BRANE1 
model is restricted to be positive. 
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The physical restrictions formulated above are necessary for a reasonable braneworld 
theory with the asymptotic limit (j3.21j) reached today, but they are certainly far from 
being sufficient for the consistency of a braneworld model with all cosmological obser- 
vations (cosmic microwave background, big-bang nucleosynthesis, etc). This important 
issue, which, in fact, constitutes the main subject of investigations in the braneworld the- 
ory, clearly lies beyond the scope of the present paper and will be considered in our future 
publications. At the moment, one can only note that the values of the constants C\ and 
C*2 in a braneworld that today evolves in regime (|3.21j) are roughly constrained by the 
values of the detected cosmological constant and gravitational coupling, respectively, as 

C^V.lHl C 2 ^ 8 -^, (3.29) 

where H pa 73 km/sec • Mpc is the current value of the Hubble parameter, and is 
Newton's constant. With these constraints, there remains a considerable freedom in the 
parameter space of the theory, which is to be further restricted by a more careful analysis. 

If condition ()3.22|) is not satisfied today, then regime ()3.21|) has not yet been reached, 
and one should turn to the general evolution equation (J2.14j) . which has a rather rich 
and interesting phenomenology for modeling dark energy For example, the current 

acceleration of the universe expansion may be described by a braneworld model which 
has C\ = in the asymptotic expansion (j3.21|) . According to this model, the current 
acceleration of the cosmological expansion is temporary, and the future evolution proceeds 
in a matter-dominated regime (see [I] for more details). Another interesting example is 
the model of loitering braneworld [23J, which can be realised even in the spatially flat 
universe. It is the negative dark radiation (the constant C in Eq. ()2.14j0 which plays 
the crucial role in the dynamics of this model and which, therefore, also needs to be 
constrained by cosmological observations. 



3.2.3 Stability of the spatially flat static solutions 

For k = 0, our equations ()3.9j) . ()3.10j) and (j3.18j) simplify to 

- 3M 6 A b + [p(l + 3w) - 2a] ■ {p + a) = , (3.30) 



and 



W 



1 + 



2m 2 
3A/ C 



(Po + o) 



-l + 3w— + \ — 
al po + a a% _ 



respectively. 

Introducing the dimensionless parameters 

p f 2m 2 p at 2 2m 2 a 



x 



6m 2 3M 6 



y 



6m 2 



3M 6 



2m 4 A b 
3M 6 



(3.31) 
(3.32) 



(3.33) 
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we can express the quantities that we used before in their terms. Since only two of these 
parameters are independent, it is convenient to express the relevant quantities in terms 
of x and y. We have 

z = ~{x + y)[{l + 3w)x-2y] , (3.34) 
C 3M 6 



[l+w)x(x + y), (3.35) 



W = ™l {1 + w)x W + *»)*- (3.36) 
4m 4 1 + x + y 

Taking into account the physical restrictions (J3.25|) and (J3.28)) and equations (J3.34j) 
and ()3.3fi|) . we can indicate the domains of existence of static solutions (both stable and 
unstable) in the (x, y) plane. 

The natural condition po > leads to the restriction x > 0. 

For BRANE1, restrictions (j3.25|) additionally give 

y>0, Ay > (x + y) [(1 + 3w)x - 2y] , (1 + 3w)x > (1 - 3w)y . (3.37) 

For BRANE2, restrictions (J3~2%|) additionally give 

y>-l, Ay + 2 > (x + y) [(1 + 3w)x - 2y] . (3.38) 

The corresponding domains are shown below in Figs. Q and 121 for the two important 
cases w = and w = 1/3. The shaded regions in these figures correspond to unstable 
static solutions (W > 0), and the dotted regions correspond to stable static solutions 
(W < 0). In general, stable static solutions exist only for the BRANE2 model in the case 
w < 1/3. 

It also follows from (J3.35|) that spatially flat static solutions with matter satisfying 
the strong energy condition can exist in the BRANE2 model with C — 0, i.e., with zero 
dark-radiation term, in contrast with the Randall-Sundrum model, where spatially flat 
or open static solutions require negative dark-radiation term. This difference is connected 
with the fact that the physically admissible range of parameters in the Randall-Sundrum 
model requires positive brane tension while, in the presence of the scalar-curvature term 
in the action for the brane, this restriction is weakened, and negative brane tensions are 
also allowed. As follows from ()3.35|) . the condition C = is equivalent to x + y — 0, 
and the region with negative y in the admissible domain of parameters allowing for this 
condition can be seen in Fig. |21 

Stability conditions for a spatially non-flat universe (k ^ 0) are further considered by 
a different method in Sec. El 
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Figure 1: BRANE1, w = (left) and w = 1/3 (right). The shaded regions indicate the physically 
allowable parameters, corresponding to unstable static solutions. 




Figure 2: BRANE2, w = (left) and w = 1/3 (right). The shaded regions of the physically 
allowable parameters correspond to unstable static solutions, and the dotted region corresponds 
to stable static solutions. 

4 Non-static solutions 

In this section, we describe non-static solutions of the braneworld theory under investiga- 
tion and also determine their stability. To do this, we use a set of convenient phase-space 
variables similar to those introduced in [Bl Ej . The critical points of the system of differ- 
ential equations in the space of these variables describe interesting non-static solutions. A 
method for evaluating the eigenvalues of the critical points of the Friedmann and Bianchi 
models was introduced by Goliath and Ellis p] and further used in the analysis by Campos 
and Sopuerta [H| of the Randall-Sundrum braneworld theory. In this section, we extend 
their investigation to the case where the scalar-curvature term is also present in the action 
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for the brane, but we restrict ourselves to the homogeneous and isotropic cosmology. 



4.1 The case of the Randall— Sundrum braneworld model 

In this subsection, we briefly describe the method and results of [H] , where the model with 
m = was under consideration. After setting m = 0, Eqs. (|2.11|) and (j2.14j) turn to 



9M b 

and 



O „ 2 , K A b 

~ + n + — - — 

a a z 6 



+ \p{l + Zw) -2a]{p + a) =0 (4.1) 



2 « _ (p + ct) 2 C A b 

respectively. Here, for simplicity, we consider the case C > 0. The last equation can be 
written as follows: 

Using the definition of the four- dimensional cosmological constant 

.A. ^ 

Ars = f + 3^ , (4-4) 

which is proportional to the Randall-Sundrum constraint [l u , we introduce the notation 
similar to those of [Bj: 

fip= 9M^' fifc = ~^' fiA= 3^' ^ A = 9mW' fic = ^' (4 ' 5) 

The authors of jH] work in the five-dimensional f2-space (Q p , Qk, &c)- However, 

the Q parameters are not all independent since they are related by the condition 

tt p + tt k + tt A + tt x + tt c = 1 . (4.6) 

Taking this constraint into account leaves us with a four- dimensional space and only four 
eigenvalues of the static solution to be found. 
Introducing the primed time derivative 

, _ 1 d 
~ Hdi : 

one obtains the system of first-order differential equations 

tt p = [2(l + q)-3(l + w)]n p , 

a A = 2(l + q)Q A , (4.7) 
tt x = 2[(l + q)-3(l+w)]Qx, 
tf c = 2(q-l)n c , 
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where 

q = -^^tt p - n A + (3w + 2)fi A + n c . (4.8) 

The behavior of this system of equations in the neighborhood of its stationary point is 
determined by the corresponding matrix of its linearization. The real parts of its eigen- 
values tell us whether the corresponding cosmological solution is stable or unstable with 
respect to the homogeneous perturbations. 

As was said above, only four parameters Q are independent because of constraint ([4.6)1 . 
In each of the cases listed below, we choose independent parameters conveniently, and our 
choice is reflected in the indices of the corresponding eigenvalues. 

(1) The Friedmann case, or (Q p , fl k , fl A , Q\, Q c ) — (1; 0, 0, 0, 0). We have 

Q = , (4-9) 

and the eigenvalues are 

X k = l + 3w, A A = 3(l + iu), A a = -3(1+u>), X c = Sw - 1 . (4.10) 

(2) The Milne case, or (Q p , Q k , Q A , Q x , Q c ) = (0, 1, 0, 0, 0). We have 

q = 0, (4.11) 

and the eigenvalues are 

A p = -(l + 3w), A A = 2, A a = -2(2 + 3w) , A c = -2 . (4.12) 

(3) The de Sitter case, or (fl p , fl k , Q\, Q c ) — (0, 0, 1, 0, 0). We have 

q=-l, (4.13) 

and the eigenvalues are 

X p = -3(l + w), A fc = -2, \ x = -6(l+w), A c = -4. (4.14) 

(4) The case considered in [21], or (fl p , fl k , fl A , Q c ) — (0, 0, 0, 1, 0). We have 

q = 2 + 3w, (4.15) 

and the eigenvalues are 

X„ = 3(l + w), X k = 2(2 + 3w), X A = Q(l + w), X c = 2(l + 3w). (4.16) 

(5) The dark-radiation case, or (fi p , fl k , fl A , Q\, Q c ) — (0, 0, 0, 0, 1). We have 

q = l, (4.17) 

and the eigenvalues are 

X p = l-3w, X k = 2, A A = 4, A a = -2(1+3w). (4.18) 
All these results for the case m = were obtained in 
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4.2 The case 

In this subsection, we apply the method of 6j to the case of braneworld theory with the 
scalar-curvature term in the action for the brane, i.e., where m ^ 0. We introduce the 
notation 



b = ^, p = ^_, s- n 



3£ 2 H 4 ' 6m 2 H 2 ' 6m 2 H 2 ' 

1 K=" D C 



(4.19) 



£ 2 H 2 ' a 2 H 2 ' £ 2 a A H A 

and express the deceleration parameter 



9 = -JrT 4 - 20 

from the second-order differential equation (J2.11j) : 

HI + K) - B + [P(l + 3m) - 2S][P + S - 1(1 + if)] 
,= L + P + S-l(l + A-) ' (4 ' 21) 

Equation ()2.13|) implies the constraint 

[1 + K - 2(P + S)f = 4L (1 + iT) - 2S - D . (4.22) 

Note that the introduced variables are infinite for H = and m = and, therefore, 
cannot describe these particular cases. The second case was already considered in the 
previous subsection and in 0, and the first case will be considered below in Sec. El (using 
different variables). 

Introducing the prime derivative 

we obtain the system of equations describing the evolution of our dimensionless parame- 

P> = -P(l + 3w -2q) , 
S' = 2S(l + q), 

L' = 2L(1 + q) , 

1 H)) (4.23) 
B' = AB(l + q), 

K' = 2Kq , 

D' = 4Dq, 

where q is given by Eq. ()4.21|) . The critical points of this system are those values of (P, 
S, L, B, K, D) which nullify its right-hand side. Since three equations in this system 
have similar form, there are three different options for the generic values of w, and there 
are also three special cases, namely, w = —1, w = 1/3 and w = —1/3, which we consider 
separately. 
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4.2.1 Generic situation 



Case 1. 1 + 3w - 2q = and S = L = B = K = D = 0. Using we obtain 

P = \- (4-24) 

Hence, 

H 2 = ^ , (4.25) 
6m z 

which describes the standard Friedmann universe and is similar to case 1 of the previous 
subsection. 

Case 2. q = and P = S = L = B = 0. 

In this case, Eq. (I4.21|) gives the identity, and the constraint equation (|4.22jl reads 

(1 + K) 2 = —D . (4.26) 

The universe described by this solution is dominated by the spatial curvature and/or dark 
radiation with negative constant C. In the absence of dark radiation, we have the spatially 
open Milne universe, similarly to case 2 of the previous subsection. 

Case 3. q = -1 and P = K = D = 0. Using P~2T| or (JP2]) . we obtain 

-B = ^ + 2S 2 -2S -2L. (4.27) 

This gives us only one equation for the three parameters L, B, and S. This condition can 
also be obtained from the first-order differential equation (J2.13)) . If we wish to specify the 
theory (BRANE1 or BRANE2) to which our solution belongs, we can use the equation 
following directly from ()2.14|) : 



l -- S - L = ±JL 2 + 2LS - l -B . (4.28) 

All solutions have H = const (the de Sitter case). Moreover, it follows from (j4.36J) 
below that all such solutions are asymptotically stable for all realistic w. 

The Hubble constant H can be obtained from the first-order equation (|2.14j) with 
C = 0, p = 0, and k = 0: 



H 2 = + 1 

3m 2 £ 2 



1± A 1 + 



a A b 



3m 2 6 

This case is similar to case 3 of the previous subsection. 



(4.29) 



We see that there are no analogues of the solution j24j (case 4 of the previous subsection) 
in the generic model with rn^O. This can be understood from the comparison of the high- 
energy limits of equations (J2.14|) and (J2.15|) . In the first case, the high-energy behaviour is 
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dominated by the usual linear term with respect to the energy density on the right-hand 
side, with negligible contribution from the part with the square root. Thus, it reproduces 
the usual Friedmannian behaviour. In the second case, corresponding to m = 0, the 
high-energy evolution is dominated by the term quadratic in the energy density. 

As for case 5 of the previous subsection describing a universe dominated by dark 
radiation, it is now replaced by case 2, in which the dark-radiation contribution evolves 
similarly to that of the spatial curvature. In the absence of dark radiation, it reduces to 
case 2 of the previous subsection, describing the Milne universe. 

4.2.2 Special cases 

Apart from the generic cases, there exist special cases w = —1 and w = —1/3, in which 
there arise an additional similarity in the form of otherwise different equations in ([4.23)1 . 

Case 1. 1 + 3w = —2, or w = —1. 

In this case, for q = — 1, P is not necessarily equal to zero. Matter of this type is 
equivalent to the four-dimensional cosmological constant (brane tension a) and can be 
effectively "eliminated" by a simple redefinition of a . 

Case 2. 1 + 3w = 0, or w = —1/3. 

This case is a little more complicated than the previous one. Now P is not necessarily 
equal to zero for q = 0. The constraint equation (|4.22j) in this case yields: 



This solution corresponds to the regime where the matter with negative pressure p = 
— p/3, the spatial-curvature term and the dark-radiation term (with negative constant C) 
evolve similarly in the cosmological equation (j2.14jl . so that one has, asymptotically as 
a — > 0, 



4.3 Stability of the generic solutions 

To study the stability of a critical point of our system (|4.23|) . we must expand the system 
in the neighborhood of this point and then calculate the eigenvalues of the corresponding 
evolution matrix. We introduce the notation: 



P = P + r, S = S + s, B = B + b, L = L + £ , K = K + k , D = D + d , 



and consider the same cases as in Sec. 14.2.11 In doing this, we note that the constraint 
equation (|4.22|) reduces the number of variables to five, and, moreover, the groups of 




(4.30) 




(4.31) 



(4.32) 
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variables {S, L, B} and {K, D} remain proportional between themselves during the 
evolution, so that the corresponding eigenvalues should also be proportional. 

Case 1. q = + „ 3 " ) , and we obtain the system 

s' = 3s(l + w), Z' = 3£{l + w), b' = 6b(l+w), k' = k(l + 3w) , d! = 2d(l + 3tw) . 

(4.33) 

Its eigenvalues are: 

As = \ L = 3(1 + \ B = 6(l + w), X K = l + 3w, X D = 2(l + 3w). (4.34) 

Case 2. q = 0. In this case, the set of stationary points is a curve in the parameter 
space; hence, one of its eigenvalues will be zero. We have the system of equations for the 
conveniently chosen variables 

r' = -r{l + 3w), s' = 2s, £' = 2£, b' = 46, k' = 2K Sq(r, s, £, b) , (4.35) 

where 5g is the linearized expression (|4.21j) for the deceleration parameter q (which does 
not depend on k). Its eigenvalues are 

A 1 = 0, \ P = -{l + 3w), \ S = \ L = 2 1 A B = 4. (4.36) 

Case 3. qo = —1. In this case, the set of stationary points is a two-dimensional 
hypersurface in the parameter space, so two of the eigenvalues will be zero. The system 
of equations is 



(4.37) 



r' = -3r(l + w) k' = -2k, d' = -4d , 

s' = 2S 8q(r, k, d, s, f ) , f ' = 2L 5q(r, k, d, s, f ) . 

Its eigenvalues are: 

A 1 = A 2 = 0, X P = -3(l + w), X K = -2, X D = -A. (4.38) 

If some of the eigenvalues are positive, then the corresponding solution is unstable. 
Summarizing the results obtained, we conclude that the Friedmann solution is a repeller 
(always unstable) for all w > — 1, and is a saddle point for —1 < w < —1/3; the dark- 
radiation static solution is a repeller for — 1 < w < —1/3, and is a saddle point for 
w > —1/3; the de Sitter solution is an attractor (always stable) for all w > — 1. 

5 Static case with nonzero spatial curvature 

To investigate this case, we define the new appropriate variables by dividing those of the 
previous section by L. Specifically, 

*-l *-§• <"> 
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Obviously, these variables are finite in the limit H — > 0. Moreover, only P, K, D and N 
vary with time while B and S are constants. We introduce also the new primed derivative 

d 



dt 



Then 



Stationarity of solution implies 



N' = £ 2 H = i 2 - - N 2 , 
a 

P = -3NP(1 + w) , 
K' = -2NK , 
D' = -AND . 

iV = 0, - = 0. 



(5.2) 



Since 



a 
a 



1 N + K-B + [p(l + 3w) -25] [P + S- l(N + X)] 



£ 2 



l + P + ,S-i(Ar + X) 



we obtain the equation 



= 0. 



(5.3) 



(5.4) 



(5.5) 



P{\ + 3w) -2S 

Introducing the perturbed variables 

N — n, K = K + k, P = P + r, D = D + d , S = S , B = B , (5.6) 
and linearizing the system with respect to n, k, d, and r, we obtain 

k(l-\ [P (l + 3w) - 25 ]) + r [2P (1 + 3iu) - 5 (1 - 3w) - ^(1 + 3w)] 



l + Po + S 



n = 

= Fk + Gr, 
k' = -2nK , 
r' — — 3nP (l + w) , 
d! = -4nP , 

or, introducing the vector notation z — (d, k,r,n), 

z' = Xz, 

where 



X = 



K 



( 











-4A, 


\ 













-2K 















-3P (l + w) 




\ 





F 


G 





/ 



(5.7) 
(5.8) 

(5.9) 
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The eigenvalue equation for this system is: 

-A 4 -2X 2 K F -3X 2 (1 + w)P Q G = 0. (5.10) 

Therefore, 

a 1j2 = o, \ 3ii = ±eVw (5.11) 

with 

2K Q [l-±(P Q (l + 3w)-2S )] 

e 2 (i + p + s -f) 

3P (1 + w) [2P (1 + 3w) - SpQ. - 3w) - f (1 + 3w)} 

P(l + Pa + S -f) ■ [bAZ) 

We note that two of the eigenvalues are equal to zero. In principle, the complicated 
expression (|5.12|) determines the regions of stability (W < 0) and instability (W > 0) of 
the static solutions of the theory. 

In the partial case k = (or K = 0) that was considered in Sec. El we obtain 

3P (1 + w) [2P (1 + 3w) + S (3w - 1)] 
W ~ ^(1 + P + S ) • (5 ' 13) 

In the limit m — > 0, our result reproduces the corresponding result of 0. Indeed, we 
have 

k£ 2 Anm A pi 2 Apm 2 Aam 2 

Then, in the limit m — > 0, 

W = ^ + ^^[2P(1 + 3") + ^(3^ - 1)] . (5.15) 

Expressing n/a 2 through p and a from (|4.1|) and (|4.2j) . 

« _ p(p + er)(l + w) 
^2 ~ 3M6 

we obtain 



(5.16) 



W = ^^[2p(2 + 3w) + a(l + 3w)] , (5.17) 



which coincides with the result of j^J. 



6 Conclusions 

In this work, we considered general cosmological solutions and their stability with respect 
to homogeneous and isotropic perturbations in the braneworld theory with the induced- 
curvature term in the action for the brane. In our approach, either the initial conditions 
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or the constants of the theory were perturbed. Part of the results are similar to those 
obtained by Campos and Sopuerta [H| for the Randall-Sundrum model. Specifically, the 
expanding de Sitter solution is an attractor, while the expanding Friedmann solution is a 
repeller in the phase space of the theory. However, in the model with m^O, there are no 
analogues of the solution considered in (case 4 of Sec. 14. 1(1 . and an expanding universe 
dominated by dark radiation (case 5 of Sec. 14. 1ft is replaced by a somewhat different regime 
(cases 2 of Sees. 14. 2. H and 14.2.2)1 . in which dark radiation contributes similarly to matter 
with negative pressure p = — p/3 and to spatial curvature. 

The possibility of static solutions, including those dominated by negative dark radia- 
tion, was previously discussed in [HI H2j in the context of the Randall-Sundrum model. In 
the braneworld theory with the scalar-curvature term in the action for the brane, static 
solutions with matter satisfying the strong energy condition exist not only with closed 
spatial geometry but also with open and flat ones even in the case where dark radiation 
is absent, as well as in the case where it is negative. This brings to attention an inter- 
esting possibility that the braneworld universe, even being spatially open or flat, could 
have passed through a quasi-static (or "loitering") phase, the details of which are further 
investigated in 
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